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ABSTRACT 

We show that binned differential luminosity functions constructed using the 
\/Va method have a significant systematic error for objects close to their parent 
sample's flux limit(s). This is particularly noticeable when luminosity functions 
are produced for a number of different redshift ranges as is common in the study 
of AGN or galaxy evolution. We present a simple method of constructing a binned 
luminosity function which overcomes this problem and has a number of other ad- 
vantages over the traditional 1/14 method. We also describe a practical method for 
comparing binned and model luminosity functions, by calculating the expectation 
values of the binned luminosity function from the model. 

Binned luminosity functions produced by the two methods are compared for 
simulated data and for the Large Bright QSO Survey (LBQS). It is shown that 
the \/Va method produces a very misleading picture of evolution in the LBQS. 
The binned luminosity function of the LBQS is then compared to a model two 
power law luminosity function undergoing pure luminosity evolution from Boyle 
et al. (1991). The comparison is made using a model luminosity function averaged 
over each redshift shell, and using the expectation values for the binned lumi- 
nosity function calculated from the model. The luminosity function averaged in 
each redshift shell gives a misleading impression that the model over predicts the 
number of QSOs at low luminosity even when model and data are consistent. The 
expectation values show that there are significant differences between model and 
data: the model overpredicts the number of low luminosity sources at both low 
and high redshift. The luminosity function does not appear to steepen relative to 
the model as redshift increases. 
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1 INTRODUCTION 

For three decades considerable effort has been spent in 
obtaining samples of AGN to investigate their cosmolog- 
ical evolution. This evolution is seen as a change in the 
luminosity function with redshift. In recent years, galaxy 
redshift surveys and X-ray cluster surveys have become 
sufficiently deep that the evolution of galaxies and clus- 
ters of galaxies can be examined in the same way as that 
of AGN. Quantifying and understanding the evolution of 
these populations is essential to modern cosmology. 

A number of methods exist for demonstrating and 
quantifying the cosmological evolution of a population, eg 
the {V/Vmax) method of Schmidt (1968) or by compari- 
son with Monte Carlo simulations (La Franca & Cristiani 
1997). The simplest, most intuitive method is to construct 
a binned differential luminosity function in a number of 
redshift intervals. For flux limited samples (which are the 



norm in evolution studies), this is usually done using the 
1/Va method (eg Maccacaro et al. 1991, Ellis et al. 1996). 
However, the 1/Va method introduces a significant error 
for objects close to the flux limit, and hence is not the op- 
timum method for constructing binned differential lumi- 
nosity functions. We present a simple, alternative method 
which we show is superior to 1/Va. 

This paper is laid out as follows. Section ^ describes 
both the \/Va method, and our new method, for construct- 
ing binned differential luminosity functions. The relative 
merits of the two methods are compared in Section ^ and 
in Section ^ a Monte Carlo simulation is used to demon- 
strate the improvement that our method offers over the 
1/Va method. In Section |5| we present a technique for com- 
paring binned and model luminosity functions. We apply 
the two different methods to construct binned luminosity 
functions of the Large Bright QSO Survey sample in Sec- 
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tion ^. In Section ^ this sample is compared to one of the 
model luminosity functions and evolution laws from Boyle 
et al. (1991) using both an averaged luminosity function 
and using the method given in Section ^ Our conclusions 
are presented in Section H. 



2 METHOD 

2.1 Fundamental quantities 

We define the differential luminosity function </!>, of any 
extragalactic population, as the number of objects per unit 
comoving volume per unit luminosity interval, i.e. 



cP{L,z) 



-{L,z) 



(1) 



dVdL' 

where A'^ is the number of objects of luminosity L found in 
comoving volume V at redshift z. We assume is a con- 
tinuous function over the range of redshift and luminosity 
for which it is defined. 

The differential luminosity function is often defined 
as a function of logarithmic luminosity or (by opti- 
cal astronomers) magnitude rather than luminosity. The 
method described here is equally applicable to luminosity 
functions defined in this way, with the appropriate substi- 
tution of 'log (L)' or 'm' where we have used 'L\ 



2.2 The 1/Va method 

The 1/Va method was originally proposed (Schmidt 1968) 
to measure space density (dN/dV); the method was gen- 
eralised for samples with multiple flux limits in Avni & 
Bahcall (1980). If A'^ objects have luminosities and red- 
shifts in the interval ALAz around the bin centre (L, z): 

N 



dN 
dV 



{L,z)^Y. 



Va{i) 



where Va{i) is the survey volume in which object i with 
luminosity L{i) could have been detected and remained 
in the redshift bin /S.z. It is an unbiased estimator of the 
space density (Felten 1976); the main advantage of this 
method (compared with simply dividing the number of 
objects found by the average volume searched) is that it 
takes account of the fact that more luminous objects are 
detectable over a larger volume than the less luminous 
objects. 

The 1/Va method is used to estimate by dividing the 
space density by the luminosity bin width. We will refer 
to the approximation of (j) obtained using this method as 



iV 



(2) 



2.3 A binned approximation to (j) 

The luminosity function as defined in equation ^ is related 
to the expected number of objects found in any region 
ALAy(A2) of the volume-luminosity plane by 



(7V> 



ML, z)—^dzdL 

dz 



(3) 



where Zmin is the bottom of the redshift interval Az and 
Zmax{L) is the highest redshift possible for an object of 
luminosity L to be detected and remain in the redshift 
interval Az. Angled brackets denote expectation value. If 
changes little over ALAz, then 



{N)^4>{L,z) 



——dzdL 

dz 



(4) 



It follows that if A'^ objects have been found over some 
volume-luminosity region A, 

N 



dzdL 



(5) 



where 4)est is our binned estimate of the luminosity func- 
tion. 



3 COMPARISON OF THE TWO METHODS 

3.1 Which is the better estimate of ? 

Consider a single bin, corresponding to the redshift - lu- 
minosity interval AL Az around the point {Li,zi), in a 
binned luminosity function constructed from a fiux limited 
sample of objects. We assume that AL and Az are suffi- 
ciently small and/or that is a sufficiently slowly varying 
function of L and z that our single estimate of is appro- 
priate over the whole interval. Some portion of the region 
AL Az may represent objects which are fainter than the 
survey flux limit. This situation is illustrated in Fig. |l| 
where the curve L = Liim{z) represents the flux limit of 
the survey. 

The shaded region in Fig. ^a is that region of the vol- 
ume luminosity plane in the interval AL Az which has 
been surveyed. This shaded area corresponds to the dou- 
ble integral in equation ^. By definition (equation ^) , (pest 
gives a good estimate to ^(Li, zi). 

The shaded region in Fig. nib has an area equal to 
(AL Va{i)) for object i, represented by the black spot. 
This area is clearly not the same as that of the surveyed 
region (the shaded region in Fig. |l^) and hence there is 
no guarantee that (pi/Va Si^e a good estimate of 0. 
However, there are two specific cases in which (pi/Va 
always give a good estimate of <j): 

1) The entire redshift, luminosity interval AL Az corre- 
sponds to objects brighter than the flux limit of the sur- 
vey, i.e. in Fig. ^ the curve L = Lii,n{z) passes above the 
AL Az region. In this case. 



dVdL = AL AV{Az) 



and 



'l/V"„ = Vest 



2) AL is very small: 

_}_dN_ 

AllV dVdL 



d^N 



as AL -> 



Note that even in these limiting cases, <l>i/Va does not give 



a better estimate to d> than 



instead 



and 



have the same value in case 1 and converge to the same 
value in case 2. Case 1 generally applies for luminosity 
functions of objects which are much brighter than the 
flux limits. Hence for luminosity functions in a number 
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Figure 1. Volume-luminosity space 'available' to an object (the black dot) in (a) (pest and (b) ipi/Va - h=hii^{z) is the 

minimum detectable luminosity of an object of redshift z and is determined by the flux limit of the survey. 



of redshift intervals, the two methods will produce the 
same results for the highest luminosity bins of any given 
redshift shell. At the lowest luminosities in each redshift 
shell, where objects are close to the survey flux limit and 
the shaded region of volume-luminosity space in Fig.^ 
is much smaller than the rectangular ALAV region, the 
two methods are likely to give the most different results. 
In all cases, 4>eat can be expected to give as good or bet- 



ter estimates for 
estimator. 



than i 



'l/Va 



and is therefore the better 



3.2 Statistical uncertainty in the two methods 



From the deflnition of (pest in Section 2.S it is easily shown 
that the statistical uncertainty, 5(j)est, is given by: 
5N 



(6) 

where 5N is the uncertainty on A'^ objects and can be cal- 
culated from Poisson or Gaussian statistics as appropriate. 



The statistical uncertainty on 



is harder to es- 



timate, because each object makes a different statistical 
contribution. The uncertainty is typically estimated by the 
following expression: 



1 

aI' 



N 



(7) 



(eg Marshall 1985, Boyle et al. 1988) but this assumes 
Gaussian statistics and therefore is not appropriate for 
bins containing few objects. 



Hence the uncertainty on (f>est is easily calculated for 
any number of objects within a bin, but the uncertainty 
on (j)i/Va can only be properly estimated when there are 
many sources per bin. 



3.3 Multiple flux limits 

The comparison so far and Fig. |^ has assumed the simple 
case of a survey with a single flux limit. However, many 
surveys have different flux limits for different parts of the 
survey area, and most luminosity function investigations 
combine different surveys with different flux limits. In- 
deed, combining large area shallow surveys and small area 
deep surveys is currently the only practical way to ob- 
tain a luminosity function which spans a wide range of 
luminosities at a wide range of redshifts. The most effi- 
cient means to combine areas with different flux limits is 
to assume that each object could be found in any of the 
survey areas for which it is brighter than the correspond- 
ing flux limit. This is 'coherent' addition of samples in 
the language of Avni & Bahcall (1980). In this case, it 
is no longer true that only the lowest luminosity objects 
in each redshift shell are close to the flux limit, because 
objects of higher luminosity are close to the flux limit of 
the shallower survey areas. The problem with 4>i/Va, there- 
fore affects luminosity bins that correspond to any of the 
survey flux limits, although the severity of the problem is 
watered down by the deeper survey areas. 

Many surveys also have a maximum flux limit beyond 
which sources are too bright to have been selected. In this 
case Zmin in equations |l| - as well as Zmax, is a func- 
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tion of L. As an illustration, if the curve L — (z) in 
Fig. |l] were to represent a maximum (rather than a min- 
imum) flux limit, objects would have to lie on the other 
side of the curve to be included in the survey. If the ob- 
ject (black dot) in Fig. |l| were shifted straight upwards 
to lie above the black line, it is now the shaded area of 
Fig. |l|a subtracted from AL AV{Az) which corresponds 
to the double integral in equation 0, and the shaded area of 
Fig. ^ subtracted from AL AV(Az). which corresponds 
to (AL Va{i))- Just as for the minimum flux limit (sec- 
tion ^.l| ), these two quantities are different, and there is 
no guarantee that (j^i/Va gi"vs ^ good estimate of (j>. 
Therefore, for surveys with both bright and faint flux lim- 



its, 



is likely to be inferior to 



at the highest as 



well as the lowest luminosity bins in each redshift shell. 



4 A SIMULATED LUMINOSITY FUNCTION 

In this section we use a sample simulated by a Monte Carlo 
method to demonstrate the improvement that (j>est offers 
over the traditional (pi/Va- The simulation was performed 
using a two power law model X-ray luminosity function 
which is unchanging with redshift (i.e. no evolution) and 
a single flux limit. The Monte Carlo simulation produced 
a source list of ~ 10000 sources. Binned luminosity func- 
tions were produced for this simulated data in a range of 
redshift intervals using both methods. These are shown in 
Fig. ^, (j>est on the left and 4'i/Va on the right. The model 
luminosity function is shown as a dashed line. 

The two methods are in exact agreement for the high 
luminosit y p oints of each redshift interval, as expected 



from Sec 3.1. However, for the lowest luminosity bins of 



the 0.2 < 2 < 0.5 and 0.5 < 2 < 1.0 redshift intervals 
4>i/Va is signiflcantly smaller than the input model, while 
0est is a good representation of the input model. For the 
lowest luminosity bin of the 2.0 < z < 2.5 redshift shell 
(jti/Va is larger than the input model. We therefore see that 
for objects close to the flux limit 4'i/Va can be too small 
or it can be too large; in simulations we have performed 
using different flux limits and luminosity functions (pi/v^ 
is more frequently too small than too large. 



any one luminosity bin. Similarly, any model luminosity 
function which is evolving with redshift can have a large 
range of values for a single luminosity over a redshift in- 
terval (although this latter effect can be overcome if the 
evolution is known or assumed a-priori, eg Mathez et al. 
1996, Kassiola & Mathez 1990). The binned luminosity 
function could be compared to the model luminosity func- 
tion evaluated at some arbitrary 'middle point' of each 
luminosity /redshift bin; alternatively some weighted aver- 
age of the model could be produced for each luminosity 
redshift bin. 

A more valid statistical approach is to compute 
the expectation value of the binned luminosity function 
from the model, for each redshift/luminosity bin. In this 
method, a single, unique, model value is produced for each 
redshift/luminosity bin, which can be compared to the cor- 
responding binned luminosity function data point. This al- 
lows the binned luminosity function to be compared with 
the model using a statistical goodness of flt test such as 
X^, and/or enables the fltting of a model to the binned 
luminosity function. This approach is illustrated in this 
section using the Monte Carlo simulated sample of sec- 
tion U. First, the formulas for expectation value are given. 



5.1 



) 



The expectation value of (jti/Va is given by: 
and L„ 



(8) 



where Lmin and Lmax are respectively the minimum and 
maximum luminosities of objects within the luminosity 
bin and P{L) is the probability density corresponding to 
an object of luminosity L. P{L) is given by 



4>{L, z)^^dz 



(9) 



"•^'^M 

where Zmin is the lower limit of the redshift shell and 
Zmax{L) is the smaller of the maximum detectable redshift 
of an object of luminosity L and the top of the redshift 
shell. (A*') is given by equation ^. Hence 



VaiL) 



<l>{L,z)^dz dL(lO) 



5 COMPARING BINNED AND MODEL 
LUMINOSITY FUNCTIONS 

The luminosity function is frequently modelled as some 
analytical function (eg a broken power law or a Schechter 
function). Models are sometimes fitted to the binned lu- 
minosity function and tested for goodness of flt using x^- 
More often the fit is made to the unbinned data by maxi- 
mum likelihood (Crawford, Jauncey & Murdoch 1970) and 
models are tested for overall goodness of fit using a one 
or two dimensional Kolmogorov Smirnov test (Press et al. 
1992). However, even in this latter case, comparison be- 
tween the binned and model luminosity function is the 
usual recourse to find out why models are rejected, i.e. at 
which redshifts and luminosities the data and model are 
inconsistent. 

In this comparison, a problem can arise that the 
model is a continuous function and (for a steep luminosity 
function) can span a rather large range of values within 



5.2 {(pest) 

From equations |and|, the expectation value of 0est 



is given by: 



<i,^dzdL 



!^dzdL 

az 



(11) 



5.3 



((/>!/ v„) and {(pest) for the simulated 
luminosity function 

In Fig. ^ {(pi/Va) and {(ftest) are shown with </>i/v_ and 
(pest of the Monte Carlo simulated sample of section Ml. Un- 
surprisingly, the binned luminosity function of the Monte 
Carlo simulated sample is consistent with the expectation 
luminosity function derived from the input model. Note 
that the systematic differences between (pi/Va and (the 
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Figure 2. Binned luminosity functions of a simulated sample of objects using (left) (pest and (right) <^j/y^. The luminosity functions 
of successive redshift shells have been shifted in the vertical direction for clarity; the input model has no evolution with redshift and 
is shown as dashed lines. 



turn down of 4>i/Va ^^w luminosities in Fig. ^) are re- 
produced by {(j>i/Va)- Hence it is in principle possible to 
compare and/or fit a model <^ to a binned (jfi/Va provided 
tfiat tfie model is first transformed into {(pi/Va,)- However, 
it is simpler to evaluate equation |l^ for arbitrary (^(L, z) 
than equation |l^, so we are again led to the conclusion 
that 0est is superior to 4'i/Va- 

6 THE LUMINOSITY FUNCTION OF THE 
LARGE BRIGHT QSO SURVEY 

To demonstrate the difference between binned luminosity 



functions produced using 



and 



, we now apply 



these two methods to real data. For this we have chosen 
the Large Bright QSO Survey (hereafter LBQS, Hewett, 
Foltz & Chaffee 1995), the current largest single sample of 
QSOs available. This is ideal because as samples become 
larger, the systematic problems with (fn/Va become larger 
than the statistical uncertainty associated with the data. 
Note that the LBQS has a maximum flux limit as well as 
minimum flux limits (see section 3.3); these are given by 
Hewett et al. (1993). 

(Iiest and <t>i/Va s-rs shown for the LBQS in Fig. ^ for 
go ~ 0.5 and Ho=50 km s~^ Mpc^^. K-correction was per- 
formed using the composite spectrum of Cristiani & Vio 
(1990). The contrast is striking: the 01/ v„ binned lumi- 
nosity function gives the misleading impression that evo- 



lution is complex and highly luminosity dependent while 
(j)est shows a much simpler picture of the evolution; the lat- 
ter could be described approximately by pure luminosity 
evolution. 



7 COMPARISON OF THE LARGE BRIGHT 
QSO SURVEY WITH A MODEL 
LUMINOSITY FUNCTION 

We now compare the binned luminosity function ^est of 
the LBQS with a model luminosity function undergoing 
pure luminosity evolution. The model chosen is the qo — 
0.5 model from Boyle et al. (1991), which is a two power 
law luminosity function described by: 

-,0.4[Mb -M6„,fc{z)]{,3+l)-,-l 



-hlO" 



'}- 

■ 6.5 X 10"'^ mag"^ Mpc"^, a = -3.9 and (5 ■ 



where ^ 

— 1.5. The model is subject to pure luminosity evolution 
such that Mt,reak{z)=-22.^ - 8.625 log il+z), untfl z = 1.9 
after which Mbreak{z) remains constant at Af(„ea/fc(1.9). 
The binned luminosity function is shown in Fig. H; the 
different redshift shells have been placed in separate panels 
for clarity. For comparison the model is shown as both a 
volume-weighted average over each redshift shell (dashed 
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Figure 3. (left): ipeat (data points) and {<f>ast) (stepped line) and (riglit): (t>\/Va (data points) {4>i/Va) (stepped line) for the Monte 
Carlo simulated sample. 



line), and as a set of expectation values {(pest) calculated 
using Equation |l^ (solid stepped line). The two arrows 
near the bottom of each plot show the position of the 
'break' luminosity of the model luminosity function at the 
top and bottom of the redshift shell. 

The advantage of comparing the binned luminosity 
function to expectation values is well demonstrated by the 
low luminosity bins of the upper four panels of Fig. ^, in 
which the continuous volume weighted luminosity function 
makes the data look far more discrepant than it really is. 

The model from Boyle et al. (1991) has already been 
compared to luminosity function of the LBQS, by Hewett 
et al. (1993). These authors used cumulative 1/Va lumi- 
nosity functions in the same redshift shells as used here; 
these were compared to the model predictions using Kol- 
mogorov Smirnov tests in each redshift shell. Their most 
important conclusion was that &t z < 2 the luminosity 
function systematically steepens relative to the model with 
increasing z, implying that pure luminosity evolution is no 
longer a viable model for QSO evolution. Similar claims 
for other surveys have been made more recently by Gold- 
schmidt & Miller (1998), La Franca & Cristiani (1997) and 
by Hawkins & Veron (1995). 

The binned luminosity function of Fig. |^ does not 
support this conclusion. Instead, significant discrepancy 
between data and model at z < 2 is always for fainter ab- 



solute magnitudes, at or below Mbreak{z), and is always in 
the sense that the model overpredicts the number of low 
luminosity sources. To quantify this, we calculated the 
separately for bins brighter, and for bins fainter Mbreak{z) 
in the four redshift shells at z < 2. These are tabulated 
in Table ^; the lowest and highest luminosity bins of the 
0.2 < z < 0.5 redshift shell and the highest luminosity bin 
of the 1.0 < 2 < 1.5 redshift shell have not been included 
because they contain too few objects for to be appro- 
priate. At absolute magnitudes brighter than Mbreak{z) 
the model is acceptable at 95%, overall and for each red- 
shift shell individually. At fainter absolute magnitudes the 
model is significantly deviant for all but the 1.0 < z < 1.5 
redshift shell, and it is notable that in this 1.0 < z < 1.5 
shell the model shows the largest deficiency relative to the 
data at high luminosities. Although the discrepancy be- 
tween data and model at low luminosities is sufficient to 
reject the model at > 7a, there is no strong evidence that 
the luminosity function steepens relative to the model as 
redshift increases. 



8 CONCLUSIONS 

We demonstrate that the 1/Va method can lead to sys- 
tematic errors when used to produce binned differential 
luminosity functions. This problem is most significant for 
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Figure 4. Binned luminosity functions of the LBQS sample: (left) cpest and (right) (pi/Va 



-30 



Table 1. Model goodness of fit at absolute magnitudes 
brighter, and similar to or fainter than, Mi,^g„j,(2) for z < 2 





Mb J < Mbreaki^i) 


Mb J > Mbreak(z) 


Redshift shell 




P* 


xV^ 


P* 


0.2 < 2 < 0.5 


2.89/2 


0.24 


58.70/3 


1.1x10-12 


0.5 < 2 < 1.0 


3.61/2 


0.16 


14.72/3 


2.1x10-3 


1.0 < 2 < 1.5 


3.72/2 


0.16 


0.68/2 


0.71 


1.5 < 2 < 2.0 


3.30/2 


0.19 


19.21/2 


6.7x10-"' 


Total 


13.52/8 


0.10 


93.31/10 


1.2x10-15 



* Probability of being higher by chance 
if the data were drawn from the model 



objects close to their parent sample's flux limit(s). As a 
result the lowest luminosity bins of 1/Va luminosity func- 
tions which arc split by redshift can be unrepresentative, 
distorting the apparent evolution of extragalactic popu- 
lations. A new method for constructing binned luminos- 
ity functions, which does not have this problem, is pre- 
sented. The improvement of this new method over the 
\/Va method is demonstrated using a Monte Carlo sim- 
ulated sample of objects. This new method also has the 
advantages that statistical uncertainty is easily estimated 
even when there are few objects per bin. 

We also present a practical method for comparing 
binned and model luminosity functions (by eye or by sta^ 



tistical test) which resolves the problem that the model lu- 
minosity function has many values within one luminosity- 
redshift bin. 

We demonstrate the difference between the 1/14 
method and the new method for constructing binned lu- 
minosity functions with the Large Bright QSO Survey 
sample. Evolution appears complex and highly luminosity 
dependent when the 1/Va method is used, but relatively 
simple in the binned luminosity function produced using 
the new method. We also use the LBQS, along with a 
model luminosity function and evolution law from Boyle 
et al. (1991), to demonstrate the advantages of our method 
for comparing model and binned luminosity functions. We 
show that the model is inconsistent with the data at low 
luminosities, but unlike Hewett et al. (1993) we do not 
find strong evidence that the luminosity function steepens 
relative to the model as redshift increases. 
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